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ABSTRACT 

The influence of a galaxy's environment on its evolution has been studied and com- 
pared extensively in the literature, although differing techniques are often used to 
define environment. Most methods fall into two broad groups: those that use nearest 
neighbours to probe the underlying density field and those that use fixed apertures. 
The differences between the two inhibit a clean comparison between analyses and leave 
open the possibility that, even with the same data, different properties are actually 
being measured. In this work we apply twenty published environment definitions to a 
common mock galaxy catalogue constrained to look like the local Universe. We find 
that nearest neighbour-based measures best probe the internal densities of high-mass 
haloes, while at low masses the inter-halo separation dominates and acts to smooth 
out local density variations. The resulting correlation also shows that nearest neigh- 
bour galaxy environment is largely independent of dark matter halo mass. Conversely, 
aperture-based methods that probe super-halo scales accurately identify high-density 
regions corresponding to high mass haloes. Both methods show how galaxies in dense 
environments tend to be redder, with the exception of the largest apertures, but these 
are the strongest at recovering the background dark matter environment. We also 
warn against using photometric redshifts to define environment in all but the densest 
regions. When considering environment there are two regimes: the 'local environment' 
internal to a halo best measured with nearest neighbour and 'large-scale environment' 
external to a halo best measured with apertures. This leads to the conclusion that 
there is no universal environment measure and the most suitable method depends on 
the scale being probed. 

Key words: methods: numerical - methods: statistical - galaxies: evolution - galax- 
ies: haloes - dark matter - large-scale structure of Universe 
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1 INTRODUCTION 

In the paradigm of hierarchical structure formation, the evo- 
lution of the primordial density field acting under gravi- 
tational instability drives dark matter to cluster and col- 
lapse into virialised objects (haloes). Such haloes provide 
the potential we lls into which baryon s fall and galaxies sub- 
sequently form (|White fc Reesiri978h . Haloes, galaxies and 
their environments also interact and merge as structure for- 
mation unfolds with time. It therefore follows that the prop- 
erties of a galaxy should be correlated with the properties 
of its host halo, and that a galaxy's environment, its host 
halo's environment, and the dark matter density field are all 
related in some measurable way. 

Such galaxy/halo/dark matter correlations with envi- 
ronment have lead to a variety of work examining the envi- 
ronmental dependence of the physics of galaxy formation, 
both theoretical and observational. Measurements of the 
galaxy two-point correlation function and halo occupation 
distribution function (HOD) have shown that more mas- 
sive, brighter, redder, and passive early-type galaxies tend 
to be more strongly clustered and hence presumably located 
in denser environments, while the reverse is true for galax- 
ies that have lower mass, are fainter, bluer and star form- 
ing (e.g.lNorberg et al.ll2002l:IZehayi et al.ll2005l:ISheth etail 



20061 ; iLi et al J 120061 ; I Tinker et al.ll2008l ; lEllison et al.ll2009l; 



Skibba & Shcth 2009; S kibba et aDl2009l : Ide la Torre et all 



2011). 



A more direct probe of the influence of environment 
is the local density field of neighbouring galaxies around 
each galaxy (defined in various ways). These techniques are 
better suited to analysing targeted halo and galaxy envi- 
ronment correlations and have proven valuable in the cur- 
rent era of large galaxy survey data sets, where galaxy 
catalogues can be simultaneously 'sliced' in multiple or- 



galaxy properties on environment (e. 


g. iKauffmann et al. 


2004; Blanton ct al. 2005i: ICroton et al.ll2005l: ICooDer et al. 


2006 ; [Baldrv et alJ 120061: IPark et ail 


120071: lElbaz et al. 


2007; Ball, Lovedav & Brunner 2008; Cowan & Ivezic 2008; 


O'Mill. Padi la & Garcia Lambas 2008 


iTasca et all 2008: 


Ellison et al. 120091). 



In undertaking any such analysis the choice of envi- 
ronmental indicator is important and no one standard has 
yet emerged. Many of the above cited papers involve dis- 
parate selection criteria, research methods and goals, mak- 
ing direct comparisons between them difficult. The defini- 
tion of environment can vary from two-point clustering and 
marked clustering statistics, to the number or luminosity 
density within a fixed spherical or cylindrical aperture, to 
the measured density enclosed by the n-th nearest neigh- 
bour. A further complication is that these methodologies 
can be performed in either two (projected) or three (red- 
shift space) dimensions. As a consequence, some analyses 
have yielded irreconcilable results. 

All methods that attempt to quantify the environment 
around a galaxy require some parameter choices. Those that 
involve a spherical or cylindrical aperture must first choose 
a fixed smoothing scale within which to measure the local 
galaxy over- or under-density. On the other hand, when en- 
vironment measures involve the n-th nearest neighbour, the 
choice of n instead becomes important. Once n is fixed, this 



statistic adapts its scale to keep the signal-to-noise constant. 
But how should one interpret a statistic that combines the 
physical processes from widely disparate scales across one 
smoothly varying curve? And how should this be compared 
with statistics that instead fix the scale along the same 
curve? 

Further complicating comparison are the selection cri- 
teria of a dataset itself, its geometry and volume, and the 
redshift and magnitude uncertainties of the galaxies in it. 
In short, the measurement of 'environment' used in vari- 
ous studies can be completely different, and environmental 
correlations should be interpreted and compared with cau- 
tion. Some environment measures can have advantages and 
disadvantages for particular research goals. A number of au- 
thors have tested and compared a few enviro n ment measures 
(e.g. [Cooper et all [20051; IWolf et~aD 120091: [Callazzi et ail 
20091; iKovac et al.ll2010l: IWilman. Zibetti fc Budavarill2010l : 



Haas. Schave fc Jeeson-Danielll2012l ). In general, while the 
environment measures are correlated, they often exhibit con- 
siderable scatter between them. 

The primary goal of this project is to compare a vari- 
ety of published environment measures using a single well 
constrained data set. For this purpose, we take a dark 
matter halo catalogue and construct a mock galaxy cata- 
logue designed to have approximately the same global sta- 
tistical propertie s as the Sloan Digital Sky Survey (SDSS; 
I York et al.ll200fjh main galaxy sample. We then are able to 
compare the galaxy environment measures to halo mass, 
dark matter density, and to each other. We also attempt 
to answer some important questions, such as: Do the differ- 
ent environment methodologies break nicely into different 
groups that optimally sample the underlying density field in 
particular ways? Do the statistics of various galaxy proper- 
ties change dramatically in different environment bins mea- 
sured in different ways? Can we find a more fundamental 
definition of environment that is measurable observation- 
ally? 

This paper is organised as follows: In Section [2] we out- 
line the mock galaxy catalogue that was generated, con- 
strained by the SDSS, and used to study environment mea- 
sures. In Section|3]we review the range of environment mea- 
sures available in the literature that are used as part of this 
study. Having established the method, Section [4] explores 
how the different measures relate to the dark matter halo 
mass, galaxy colour and large-scale dark matter environment 
for each galaxy. We also explore how the measures relate to 
each other for an individual galaxy. Finally in Section [5] we 
discuss and summarise our findings. 



2 GALAXY AND HALO CATALOGUES 
2.1 The Millennium Dark Matter Simulation 



We b egin with the Millennium Simulation (jSpringel et al.l 
2005) which is a large A^-body simulation of dark matter 
structure in a cosmological volume. The Millennium S im- 
ulation uses the GADGET Tree-PM code l|Springelll2005h to 
trace the evolution of 10 billion dark matter particles across 
cosmic time in a cubic box of 500/i _1 Mpc on a side, with a 
halo mass resolution of ~ 5 x 10 10 7i _1 Mq. It adopts the con- 
cordance ACDM cosmological parameters, chosen to agree 
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Figure 1. The mean number of galaxies above a given luminosity 
present in dark matter haloes of different mass. Error bars denote 
the 16th and 84th percentiles and are plotted for haloes that 
on average host at least one galaxy. Lines represent the input 
model and correspond to Equation [T] [Note this Figure differs 
from the published version as this one contains all haloes, while 
the published version contains only haloes hosting at least one 
galaxy.] 



with a combined analysis of the Two-Degree Field Galaxy 
Redshift Survey r2dFGRS; IColless et alJbOQlD and the first- 
year Wilkinson Mic rowave Anisotropy Probe data (WMAP; 
ISpergel et ai]|2003l V. Q = 0.25, fi A = 0.75, h = 0.73, n = l, 
and <T8 = 0.9. 

The haloes are found by a two-step procedure. First, 
all collapsed haloes with at least 20 particles are identified 
using a standard friends-of-friends group-finder with link- 
ing parameter b — 0.2. Then, post-processing with the sub- 
structure algorithm SUBFIND (|Springel et alj 1200 it ) subdi- 
vides each friends-of-friends halo into a set of self-bound sub- 
haloes. We note that comparable halo properties ar e found 
using other structure finders fsee lKnebe et alj|201lf ). 



2.2 Embedding Galaxies in Haloes 

From the Millennium Simulation halo merger tree at z — 0, 
we construct a mock gala xy catalogue using t he halo occupa- 
tion method_d«scnbed in Skibba^eFal] (|200g, hereafter S06) 
and ISkibba fc Shetbl (|2009l . hereafter SS09); we refer the 
reader to these papers for details. Other halo-model descrip- 
tions ofj£alaxy_£fu^t^^ luminosity functions 



(e.g. lYang, Mo fc van den Bosch] [2 003) and subhalo abun 



dance matching (e.g. Kravtsov et al.l 20041 ) — would produce 



similar mock catalogues, although an advantage of the SS09 
approach is that it includes a strongly constrained model 
of galaxy colours. S06 describes how the luminosities and 
real-space and redshift-space galaxy positions are modelled. 

Our model distinguishes between the 'central' galaxy in 
a halo and all the other galaxies ('satellites'). We assume 



that central galaxies have the same positions and velocities 
as the haloes in the dark matter simulation. In other words, 
central galaxies are at the centre of the haloes, and the satel- 
lites are located around them. An important assumption 
in the model is that all galaxy properties — their numbers, 
spatial distributions, velocities, luminosities, and colours — 
are determined by halo mass alone. These galaxy properties 
are constrained by SPS S observatio ns, including the lumi- 
nosity function dBlanton et al. 20031) . luminosity-dependent 
two-point clustering ( Zehayi et al.ll2005l ; ISkibba et al"]|2006l ; 



IZheng. Coil fc Z chavi 2007), and the colour- magnitude dis - 
tribution and colour-dependent clustering IjSkibbal 12009). 
Note that the clustering constraints result in a mock cata- 
logue that approximately reproduces the observed environ- 
mental dependence of luminosity and colour, on scales of 
100 /i _1 kpc to 30 /i _1 Mpc. 

The number of satellite galaxies in the model follows 
a Poisson distribution with a mean value that increases 
with halo mass. The satellites ar e distributed around the 
halo c entre so that they follow a iNavarro. Frenk fc White] 
Jl99fj) profile with the mass-con ce ntratio n relation from 
iMaccio. Dutton fc van den Boschl ||2008h . We assign 
redshift-space coordinates to the mock galaxies assuming 
that a galaxy's velocity is given by the sum of the velocity 
of its parent halo plus a virial motion contribution that 
is drawn from a Maxwell-Boltzmann distribution with 
dispersion that depends on halo mass (S06). 

We specify a minimum r-band luminosity for the galax- 
ies in the catalogue, M r — 51og(ft) = —19, to stay well above 
the resolution limit of the Millennium Simulation, avoiding 
any issues of completeness that may bias our results. We 
generate luminosities for the central galaxies, while account- 
ing for the stochasticity between their luminosities and host 
halo mass, and then we generate the satellite luminosities so 
that the observed luminosity distribution is reproduced for 
M r -51og(ft) < -19 (S06). 

We model the observed g — r colour distribution at a 
given luminosity as the sum of two Gaussian components, 
commonly referred to as the 'blue cloud' and 'red sequence'. 
Our colour model has five constraints as a function of lu- 
minosity: the mean and scatter of the red sequence, mean 
and scatter of the blue cloud, and the blue fraction. We as- 
sume that the colour distribution at fixed luminosity is ap- 
proximately independent of halo mass, and that the satellite 
colour distribution varies such that its mean increases with 
luminosity (i.e., the satellite red fraction increases with lu- 
minosity in a particular way). These two assumptio ns are 
tested and verified with galaxy group catalogues in ISkibbal 
(Hoi). 

This procedure produces a mock galaxy catalogue con- 
taining 1.84 million galaxies, of which 29 percent are satel- 
lites. Galaxies occupy haloes with masses ranging from 10 11 
to 10 15 ' 3 /i _1 Mq. We also construct a mock light cone from 
the catalogue by selecting galaxies that are within a radial 
distance of 500/i _1 Mpc from one corner of the box. This 
gives an opening angle of 90 x 90 degrees and a depth of 
500/i _1 Mpc, for which right ascension and declinations are 
determined. The analysis in Section [3] is carried out using 
a sample of galaxies that are common to both the box and 
the cone and are chosen so not to be affected by edges. Fig- 
ure [T] shows the mean number of galaxies as a function of 
halo mass, for two luminosity thresholds (L m i n ). By con- 
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Figure 2. Left: the r-band luminosity function for the mock galaxy ca t alogue created using the HOD of ISkibba fc Sheth! j200Sf ) 
(red line) compared with that o f the semi-analytic |Pe Lucia fc Blaizod (|2007f ) model (blue line) and the SDSS observed values 
(Yan g, Mo fc van den Bosch 2009) (black points with errors). Right: Two-point correlatio n function of all, re d, and blue galaxies in 
the mock catalogue (lines), compared with the equivalent observed results in the SDSS from lZehavi et aP l|2005l ) (points with errors). 



struction, the number of galaxies consists of the number of 
central galaxies plus the number of satellites, such that 



(iV ga i|A/,i min ) = (N ccn \M,L min ) 
where, 



l + (N Bat \M, L min > 



(iV cen |M) = - 



1 + erf 



log(M/M min ) 

ClogM 



and 



M - M 



(1) 



(2) 



(3) 



(See Appendix A2 of SS09 for details). All of the free pa- 
rameters depend on luminosity. The slope of the power law, 
a, is nearly unity. One may define a parameter Mi, which 
is equ al to or slightly larger than M 1 (1 Zheng. Coil fc Zehavil 
2007), and is proportional to the minimum halo mass: Mi ss 
20M m i n . This determines the mass above which haloes typ- 
ically host at least one satellite galaxy. Therefore, since for 
M r ^ —19 the minimum halo mass is w 10 11 ' 5 h -1 M©, 
the mean number of galaxies rises rapidly like a linear 
power law at masses larger than twenty times this value, 
or « W 12 - 8 h~ 1 M Q , as seen in Figured] 

At the high halo mass end, galaxy number shows a near 
linear relationship with dark matter halo mass, which occurs 
by construction in the halo occupation model. This implies 
that the number of galaxies per unit dark matter mass is 
constant, or put another way, each galaxy contributes the 
same mass of dark matt er to the cluste r . Thi s is in agree- 
ment with the findings of lPoggianti et al l (|2010l ) and to some 



degree is the natural consequenc e of a structure bui l t hier - 
archically. This also agrees with lBlanton~fc Berlindl (|2007T) 
who find that galaxy distributions are only affected by the 
host dark matter halo, and not by the surrounding density 
field, for the SDSS galaxy group catalogue. 

The r— band luminosity function of galaxies in the mock 
catalogue is shown in the left-hand panel of Figure [5] and 
is co mpared to both the observed SD SS luminosity func- 
tion (|Yang. Mo fc van den Bos ch 2009) and a popular sem i- 
analytic galaxy formation model (IDe Lucia fc Bla izot 2007). 
In the right-hand panel of Figure [2] we show the mock two- 
point correlation functions of all, red, and blue galaxies and 
compare them with the equivalent SDSS measurements of 
IZehavi et al.l (|2005). Note t hat the colou r -depe ndent two- 
point function measured bv lZehavi et al.l (|201ll ) is slightly 
different from that constrained in the mock, likely due to 
the presence of the Sloan Great Wall in the real data, an 
unusually massive supercluster at z ~ 0.08. 

We have made the mock galaxy catalogue as real- 
istic as possible, and although the catalogue reproduces 
the observed environmental dependence of luminosity and 
colour, there are nonetheless a few limitations to the model. 
For example, we have assumed virialised (dynamically re- 
laxed) dark matter haloes even though some haloes are 
not, such as those hav ing recently experienced a merger 
(e.g. lMaccio et al.ll2007T ). We have also assumed that central 
galaxies are always the brightest galaxy in a halo and lie 
at the centre of their potential well, although in a nonzero 
fraction of haloe s, especially massi ve haloes, this assump- 
tion is not valid |Skibba et al.|[20Tll ) . Finally, we force satel- 
lite galaxy properties to depend only on halo mass, not on 
halo-centric position, although there is evidence of such a 
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dependence at fixe d mass (e.g. Ivan den Bosch et al. 2008a; 
lHansen et al.l [20091 ). While our mock galaxy catalogue re- 
sembles a spectroscopic catalogue, some environment mea- 
sure s used in the literatu re are based on photometric data 
(e.g. iGallazzi et aH [2009); for tests with such measures one 
can add scatter to the redshifted mock galaxy positions, for 
example. 



3 ENVIRONMENTAL MEASURES 

There are many different methods of measuring galaxy envi- 
ronment available in the literature. Most of these can be cat- 
egorised into two broad groups: those which use neighbour 
finding and those that use a fixed aperture. An overview of 
the methods used in this work are presented in the follow- 
ing subsections and summarised in Table \T\ along with the 
authors who implemented them. 



3.1 Nearest Neighbour Environment Measures 

The principle of nearest neighbour finding is that galaxies 
with closer neighbours are in denser environments. To cre- 
ate a standard measure for this, a value of n is chosen that 
specifies the number of neighbours around the point of in- 
terest. In its simplest form, the projected surface density of 
galaxies, <r n , can then be defined as 

o n = — , (4) 

where n is the number of neighbours within the projected 
distance r n , the radius to the n-th nearest neighbour. One 
disadvantage of quantifying environment using projected 
statistics is that two galaxies can appear close together when 
they are in fact just a chance alignment and are actually 
separated by a larger distance in the third dimension. While 
there is no simple way to overcome this observationally, one 
can adopt a velocity cut about each galaxy, typically of order 
±1000 km s" 1 , to minimise the number of such alignments. 

For data where a third dimension has been measured for 
each galaxy (e.g. redshift), the denominator of Equation [3] 
is replaced by the enclosed volume: 



= C- 



(4/3)7T7i 



(5) 



When using three dimensions careful consideration of red- 
shift distortions are needed and this often leads to two di- 
mensional projected distances often being used. The nearest 
neighbour estimator was recently ap plied to the Galaxy and 
Mass Assemb l y cat alogue (GAMA; iDriver et all l201ll ) by 
iBrough et all feoilb using the distance to the first nearest 
neighbour above a given luminosity, although typically 3-10 
neighbours are used. 

Variations on the n-th nearest neighbour approach have 
been proposed in an attempt to improve the robustness of 
statistic as a measure of l ocal density. One such method 
used bv lBaldrv et~ai] (|2006n was to take the average of two 
different neighbour densities, in their case the 4th and 5th 
nearest neig hbour projected surface densities. An alternative 
proposed bv lCowan fc Ivezicl (|2008T ) was to use the distance 
to every neighbour up to the tenth instead of just the dis- 
tance to the tenth to calculate the density. They adopted a 
Bayesian metric such that 



(6) 



where C = 11.48 is empirically determined so that the mean 
of <f> matches the number density when the density is esti- 
mated on a regular grid for a uniform field, and di is the 
distance to neighbour i. 

One can also use numerical simulations to guide the 
nearest neighbour calibration. Calculating densities using 
neighbours has long been used in Smooth Particle Hydrody- 
namics (SPH) and this technique can be applied to galaxies 
in simulations. SPH calculates the density around a point 
by weighting each neighbour based on its distance from the 
point, with the smoothed galaxy density defined as 



P = J2w(\n\,h) 



(7) 



Here, n is the number of neighbours used and V^dr^l, ft) is 
the weighting given by 



W(r, h) 



r/i :i 



1 - ( 

2(1 





< - < - 

U h 2 

- < - < 1 

2 h 5=5 1 

T > 1 



(8) 



where r is the distance to each neighbour and h is the dis- 
tance to the n-th nearest ne ighbour. This weighting corre- 
sponds to the spline kernel of lMonaghan fc Lattanziq (|l985T ) 
and is the standard kernel of SPFr. Thi s method was used 

3jT l2007l 



with 20 neighbours in iPark et all (j2007T ). but values of 32 
and 64 are more common in SPH. 

Another way to const rain loca l galax y density using 
neighbours was proposed bv lLi et "all (|201 lh f or the Redshift 
One L DSS-3 Emissio n line Survey (ROLES: iGilbank et all 
2010). Li et alj (|201 lh considered the volume element of the 
nearest neighbour found by constructing a three dimensional 
cylinder using the five nearest neighbours to define its radius 
and depth. In other words, this technique encloses the five 
nearest neighbours in a cylinder that no longer has to be 
centred on the galaxy being sampled, and leads to a better 
estimate of the relevant volume when compared with simply 
using a sphere of radius the fifth nearest neighbour. 

Further consideration of the volume can be made by 
calculating the Voronoi volumes ar ound each galaxy as 
a measure of the environment (e.g. iMarinoni et al.1 1200*21 : 
ICooper et ai1l2005lh Voronoi volumes are polyhedrons con- 
structed by bisecting the distance vectors to the nearest 
neighbours. Each galaxy will have a volume around it, for 
which it does not have to be at the centre, defining the 
points in space that are closer to it than any other galaxy. 
This gives an estimate of the local density. Unlike the other 
neighbour-based methods, the number of neighbours used 
to define the shape of the volume probed is not fixed, which 
makes the technique fully adaptive. For this study a pro- 
jected Voronoi measurement is made by collapsing galax- 
ies into two dimensional slices of 50/i _1 Mpc in depth. The 
Voronoi shapes are then constructed on these surfaces to 
calculate the surface density of each galaxy. 



1 We have adopted the notation of h corresponding to the point 
at which the kernel equals zero as opposed to 2h as is used in 
traditional SPH literature. This is just a notational change. 
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Num. 



Method 



Author 



Neighbours 

1 3rd Nearest Neighbour 

2 Projected Voronoi 

3 Mean 4th & 5th Nearest Neighbour 

4 5 Neighbour Cylinder 

5 7th Projected Nearest Neighbour 

6 10 Neighbour Bayesian Metric 

7 20 Neighbour Smooth Density 

8 64 Neighbour Smooth Density 



Muldrew 
Podgorzec & Gray 
Baldry 1 
Li 2 
Ann 
Cowan 3 
Choi & Park 4 
Pcarce 



Aperture 

9 lft-'Mpc (±1000 km s" 1 ) 

10 2h" 1 Mpc (iSOOkms" 1 ) 

11 2h~ 1 Mpc (ilOOOkms- 1 ) 

12 2h~ 1 Mpc (ieOOOkms^ 1 ) 

13 Sh-iMpc (ilOOOkms" 1 ) 

14 8 h _1 Mpc Spherical 



Griitzbauch & Conselice 5 

Gallazzi 6 
Griitzbauch &; Conselice 

Gallazzi 6 
Griitzbauch & Conselice 

Croton 7 



15 
16 
17 
18 
19 
20 



Annulus 

0.5 - 1.0/i _1 Mpc (±1000 kms" 1 ) 
0.5 - 2.0h" 1 Mpc (±1000 kms" 1 ) 
0.5 - 3.0/i _1 Mpc (±1000 kms" 1 ) 
1.0 - 2.0h- 1 Mpc (±1000 kms" 1 ) 
1.0-3.0/i _1 Mpc (±1000 kms" 1 ) 
2.0 - 3.0/i _1 Mpc (±1000 kms" 1 ) 



Wilman 
Wilman 
Wilman 
Wilman 
Wilman 
Wilman 



& Zibetti 8 
& Zibetti 8 
& Zibetti 8 
& Zibetti 8 
& Zibetti 8 
& Zibetti 8 



Table 1. List of environment measures used in this study and the a uth ors who implem ente d them, including refer ence s where applicable . 
S ee Section[3]for fur t her d etail s. References: l :lBaldrv et al.l 12006]), 2:lLi et"afl fcOllT), 3:ICowan fc Ivezid feOOSl). 4:|Park et all 12007ft . 



5: lGriitzbauch et al.l l l201ll) , 6: iGallazzi et al. I J2009h . 7: ICroton et al. I <2005ft and 8: lwiIman"*zTbetti^r^udavaril 



3 Neighbours 7 Proj. Neighbours 10 Neigh. Bayesian 64 Neigh. Smooth 




log(M H alo/M Q ) 

Figure 3. The abundance of galaxies that have environments of a given percentage rank plotted against host halo mass, where envi- 
ronment is defined by the (from left to right) 3rd nearest neighbour, 7th projected nearest neighbour, 10 neighbour Bayesian metric 
and 64 neighbour kernel smoothed (SPH style). Contours are linearly spaced showing regions of constant galaxy number. The bimodal 
distribution is caused by the neighbour search remaining in or leaving the halo to find the next nearest neighbour. 



In Section 0] we apply a number of the above nearest 
neighbour methods to the mock galaxy catalogue described 
in Section[2]and quantify their relative strengths, weaknesses 
and optimal applications. 

3.2 Fixed Aperture Environment Measures 

In contrast to nearest neighbour methods, which define envi- 
ronment using a varying scale around each galaxy set by the 
distance to a pre-determined number of galaxy neighbours, 
fixed aperture methods instead probe a fixed area or volume 
around each galaxy, within which the number of neighbours 



are counted. The more galaxies inside this area or volume, 
the denser the environment is assumed to be, and vice versa. 

Fixed aperture measures are often expressed as a den- 
sity contrast, 5, instead of a density, p. Density contrast 
rescales the aperture count with respect to the mean and is 
typically defined as 

where N g is the number of galaxies found in the aperture, 
and N g is the mean number of galaxies that would be ex- 
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pected in the aperture if galaxies were instead distributed 

randomly throughout the entire volume. 

The fixed aperture technique was used bv lCroton et al.l 

(|2005l ) to investigate the environme nts around galaxie s 
in the 2dF Galaxy Redshift Survey (|Colless et al] l2003h . 
ICroton et all (|2005l ) used spherical apertures of radius 
8ft _1 Mpc, having investigated a range of sizes from of 
4/i _1 Mpc to 12/i _1 Mpc (see also I Abbas fc Shetbll2006l ) . 

When distance information is not of sufficient accuracy 
(or absent), apertures in this methodology are instead pro- 
jected on to the sky. Where possible, authors will then im- 
pose a velocit y cut of order ±1000 km s _1 to minimise in- 
terlopers (e.g. lGriitzbauch et alJl201lT ). for the same reasons 
discussed in Section 13.11 The magnitude of this velocity cut 
can vary dep ending on distance un certainties. This was in- 
vestigated bv lGallazzi et all (|2009T > who found velocity cuts 
of ±6000 km s _1 (dz — 0.02) represent the typical photomet- 
ric redshift uncertainty and ±500 km s _1 (dz — 0.0015) rep- 
resent the typical spectroscopic redshift uncertainty. Such 
errors can often have a detrimental effect on the measured 
density if not appropriately accounted for. Note that when 
a velocity cut is imposed, an otherwise spherical aperture 
elongates into a cylinder in three dimensional space, within 
which galaxy counts are then taken. Whether this distortion 
is important for the environment measure depends on the 
focus of the analysis. Typical scales for the radius of an aper- 
ture range from lft _1 Mpc to 10/i -1 Mpc, probing environ- 
ments spanning individual haloes to large super-structures 
and voids in the cosmic web. 

A variation on the fixed ape r ture m ethod was proposed 
in IWilman. Zibetti fc Budavaril (|2010| ) , where counts were 
taken in annuli of increasing inner and outer radius, rather 
than within a single fixed aperture volume. This technique 
enables the larger scale environment to be probed and the 
influence of local regions around individual galaxies to be 
removed. In its optimal form different sized annuli are ap- 
plied in combination with apertures to better constrain the 
halo size and changes of environment with distance from the 
galaxy. 

Finally, in addition to environment being defined by 
galaxy positions within the volume, we also measure envi- 
ronment as inferred from the background dark matter distri- 
bution. To obtain the neighbourhood dark matter environ- 
ment in the Millennium Simulation the full volume is broken 
into a three dimensional grid with side- length 2 h~ Mpc. At 
the centre of each grid element a three dimensional Gaussian 
density is calculated using the local dark matter particles, 
smoothed over three different scales: 2.5, 5, and 10/i _1 Mpc. 
This Gaussian smoothed density is similar to the kernel 
smoothed density of Equation [JJ but with a dark matter 
particle mass term in the sum. 

In Section[4]we apply a number of fixed aperture meth- 
ods to the mock galaxy catalogue and measure local density 
around each galaxy. This allows us to quantify the properties 
that aperture measured densities best probe, and compare 
with the previously described nearest neighbour estimators. 



4 RESULTS 

To investigate the different properties of each galaxy en- 
vironment measure, in this section we consider how they 




15.5 



log(M Ha io/M Q ) 

Figure 4. The percentage rank of galaxy environments plotted 
against dark matter halo mass, as in Figure [3] this time for the 
Voronoi method. Contours are linearly spaced showing regions 
of constant galaxy number. Vertical lines represent typical dark 
matter halo masses that host 1, 10 and 100 galaxies with M r — 
51og(h) ^ —19 (sec Figurc[TJ 

correlate with (1) the host dark matter halo mass, (2) the 
underlying dark matter environment, and (3) the colour of 
the galaxies. 

To facilitate this we have converted the output of each 
to a 'percentage rank' for each galaxy. This is computed 
by listing the galaxies in order of increasing density, then 
assigning them a percentage based on where they appear 
in that list, with zero percent being the least dense and 
one hundred percent the most dense. Therefore, a galaxy 
with a percentage rank of ninety five has five percent of the 
galaxies in the sample denser than it and ninety five percent 
less dense than it. This normalisation provides a fairer com- 
parison between environment estimators and probes their 
relative rather than absolute distributions across the envi- 
ronment spectrum, which would otherwise be definition de- 
pendent. 

4.1 Dark Matter Halo Mass 

By design, the most fundamental property for a galaxy 
within our model is its dark matter halo mass. Halo mass 
determines both the spatial distribution of the galaxy popu- 
lation and the individual galaxy properties. Therefore, each 
environment measure should reveal some underlying corre- 
lation. Typically halo masses of ~ 10 12 h' 1 Mq correspond 
to the field, ~ 10 13 5 h^M® to groups and ~ 10 15 /i _1 M© 
to clusters. 

4-1.1 Nearest neighbour results 

Figure [3] shows contours of the abundance of galaxies 
that have environments of a given percentage rank plot- 
ted against the host halo mass, for four different nearest 
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Figure 5. The percentage rank of galaxy environment against dark matter halo mass, as in Figure(3] for (from left to right) a 1 7i — 1 Mpc 
aperture with a velocity cut of ±1000kms~ 1 , a 2/i _1 Mpc aperture with a velocity cut of ±6000 km s — 1 , a 2/i — 1 Mpc aperture with a 
velocity cut of ±500 kms -1 and a 5 h~ 1 WLpc aperture with a velocity cut of ±1000 km s — 1 . Contours are linearly spaced showing regions 
of constant galaxy number. 



neighbour-based techniques, with the number of neighbours 
increasing from left to right. These are: the 3rd nearest 
neighbour density in three dimensions, the surface density 
for the projected 7th nearest neighbour, the three dimen- 
sional density using a 10 neighbour Bayesian metric, and 
the smooth kernel three dimensional density using 64 neigh- 
bours. 

The most noticeable feature of all panels in Figure [3] 
is that galaxies divide into two distinct groups, with the 
top ~ 20 percent dense environments occupied by galaxies 
in haloes more massive than ~ 10 12 ' 5 K~ 1 Mq, and the re- 
maining ~ 80 percent of environments occupied by galaxies 
in haloes with masses lower than ~ 10 12 ' 5 fc _1 Mg. This bi- 
modality arises from the assumed association between galax- 
ies and dark matter haloes required to fit the observed lumi- 
nosity function and clustering observations, and is explored 
further below. 

Looking in more detail, the lower 80 percent of rank- 
ordered densities in Figure [3] shows no trend with halo mass, 
and as such, the term 'local environment' no longer applies. 
In terms of a characteristic halo mass for a given environ- 
ment, this result leaves individual galaxies near clusters in- 
distinguishable from isolated galaxies in voids. 

In contrast, the behaviour of the high density-halo mass 
correlation depends on the neighbour method employed. 
In the highest 20 percent environments, low n neighbour 
searches smooth away any density dependence with halo 
mass. This can be seen by comparing the far left panel 
in Figure [3] (low n) with the far right panel (high n). As 
the number of neighbours used to define environment is in- 
creased, galaxies belonging to increasingly massive haloes 
(which host an increasing number of satellites) will be la- 
belled as increasingly dense. Thus, to more precisely draw 
out the high density-halo mass environment correlations us- 
ing nearest neighbour methods, a high n is desirable. 

The first two panels of Figure [3] provide an additional 
test of the importance of projection effects. Here, the 3rd 
nearest neighbour count is performed using three dimen- 
sional redshift space distances while the 7th nearest neigh- 
bour is performed with projected galaxy positions on the 
two dimensional sky. Both methods show the same overall 
trend with halo mass. We find that, in general, projecting 



the galaxy positions simply blurs the edges of the two clouds 
with the overall shape preserved. 

Another popular neighbour-based method used for mea- 
suring environment is Voronoi volumes, as discussed in Sec- 
tion (3T] Figure [4] shows how a Voronoi defined environment 
estimator also correlates with dark matter halo mass. We see 
a similar trend to that of the other neighbour-based meth- 
ods, with the overall result close to the 7th nearest projected 
neighbour method shown in the second panel of Figure [3] 

A comparison of Figure [T] with Figures [3] and [4] reveals 
the origin of the bi-modality. Galaxies identified to be in the 
upper 20 percent dense environments tend to be those whose 
neighbour search stays within the dark matter halo due to 
a large satellite population. Such haloes are almost always 
more massive than ~ 10 12 ' 5 h _1 M@. In contrast, the lower 
80 percent density environments are identified by neigh- 
bour searches that extend beyond the halo due to a low or 
zero satellite population of significance. In general, haloes 
with few satellites almost always have masses smaller than 
~ 10 12 ' 5 /i _1 Mq, and neighbour searches will then tend to 
probe the inter-halo rather than inter-galaxy separations. 

4-1-2 Fixed aperture results 

Many authors have employed fixed apertures to probe the 
local density around galaxies, as described in Section ^, 21 In 
a similar vein to Figure [3] Figure [5] shows how various aper- 
ture sizes correlate with host dark matter halo mass when 
a projected fixed aperture is employed with a cut in veloc- 
ity space around each galaxy. In addition, the central two 
panels show how the density-halo mass correlation changes 
if the velocity cut is increased for the same sized aperture. 
This roughly corresponds to the difference one would expect 
with data h aving photometric vs. spectroscopic redshifts, as 
discussed in lGallazzi et al] (|2009l ). 

The projected fixed aperture technique yields both sim- 
ilar and different trends when compared with the near- 
est neighbour technique shown in Figure [3] The overall 
shape is the same, with galaxies in haloes of mass less than 
~ 10 12 ' 5 1i~ 1 Mq showing little correlation of halo mass with 
environment. At the high mass end there is a plume of in- 
creasing density that is much better defined than found with 
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the nearest neighbour method (especially when compared 
to choices of low n). This suggests that the fixed aperture 
methodology is a better probe of halo mass, especially for 
small apertures and velocity cuts. There is however contam- 
ination at a fixed density from low mass haloes due to their 
close proximity to the high mass halo. 

In particular, when there are enough galaxies to de- 
fine the local large-scale structure, a fixed-scale environment 
probe is much more sensitive to the power-law nature of the 
two-point correlation function, where the abundance of close 
pairs falls off rapidly beyond the halo radius. This leads to 
the galaxy count in the fixed aperture also falling off rapidly. 
In contrast, nearest neighbour environment methods adapt 
the scale probed to keep signal-to-noise fixed. Hence, the di- 
vision between a halo's interior and exterior becomes much 
less prominent. 

At intermediate to low masses there is no relation be- 
tween fixed aperture measured density and halo mass, and 
so the environment parameter breaks down, as is also the 
case for nearest neighbour environment parameters. From an 
environment point-of-view, such haloes, which usually host 
galaxy groups, may be difficult to distinguish from cluster 
outskirts and from unassociated lower-mass haloes. 

As the aperture is increased in size, the trend with halo 
mass fades when the aperture becomes much larger than 
the structures present. For example, a super-cluster with 
a collective mass of I0 16 /i _1 M would have a radiufl of 
~ 3.5 /i _1 Mpc, smaller than the 5/i -1 Mpc aperture shown 
in the far right panel of Figure[S] When an aperture becomes 
large enough the contribution of individual haloes and struc- 
tures blur and the environment-halo mass trend weakens or 
disappears. Hence, aperture size should be chosen carefully 
from the outset and be appropriate for the science questions 
of interest. 

Finally, the two central panels of Figure [S] illustrate the 
importance of velocity (or equivalently distance) uncertain- 
ties on the environment measure. Large velocity cuts, as is 
typically required with photometric data, make measuring 
environment with a fixed aperture ineffective. This occurs 
for the same reason as using large apertures. There, the 
aperture was wider than the structures of interest which 
smoothed out the signal, while here, the depth of the aper- 
ture scatters in superfluous counts from foreground and 
background objects, diluting any correlation. This does not 
apply to the highest mass clusters as they dominate the 
depth reducing the effect of interlopers. Furthermore, any 
use of the angular correlation function as a probe of envi- 
ronment must first consider the redshift distribution of the 
galaxies and the uncertainties must be well understood (e.g. 
ICoil et al.ll2003 ; iQuadri et ai]|200St) . 



4.2 Galaxy Colour 

Galaxy colour has been shown to correlate with local galaxy 
density, with galaxies in over-dense environments being red- 
der compared with tho se in under-dense environments (cf. 
clust e r and field) (e.g. Lewis et alj 20021 ; Kauffmann et all 



12004 ICooper et atll2006l ; iGallazzi et al J 120091 ) . The model 



2 Radius here is determined by finding the scale at which the 
enclosed density is 200 times the critical density of the Universe. 



we employ in this paper has a constrained global g — r colour 
distribution that mim ics that of local galaxies in the SDSS 
|Skibbafc Shetrj2009| ). Hence, the degree to which different 
environment metrics can recover this relation can be tested. 

Figure [6] shows histograms of the g — r colour distribu- 
tion for the 20 percent most dense and 20 percent least dense 
galaxies defined with the same four nearest neighbour meth- 
ods used in Figure [3] the 3rd nearest neighbour density in 
three dimensions, the projected 7th nearest neighbour, den- 
sity defined from a 10 neighbour Bayesian metric, and the 
smooth kernel density using 64 neighbours. In the 20 per- 
cent most dense environments, all nearest neighbour-based 
environment measures show a clear red peak and a more 
weakly populated blue cloud. In contrast, in the lowest 20 
percent of environments galaxies are split more evenly be- 
tween the red and blue populations. As the neighbour num- 
ber is increased (from left to right), there are only small 
changes in the relative colour distributions in environment 
extremes. The Kolmogorov-Smirnov probability that both 
samples are drawn from the same distribution is shown in 
the upper right of each panel. 

Figure [7] shows histograms of colour for the 20 percent 
most dense and 20 percent least dense galaxies as probed 
by fixed apertures of various size, as used previously in Fig- 
ure [S] The central two panels show how these distributions 
change if the velocity cut is increased or decreased for the 
same sized aperture. This roughly corresponds to the dif- 
ference bet ween photometric an d spectroscopic redshift un- 
certainties (|Gallazzi et alj|2009h (see Section |33J|. For small 
apertures, the colour distributions of both density extremes 
look remarkably similar to that found for the nearest neigh- 
bour methods shown in Figure [6] However, as the volume of 
the fixed aperture is increased similar trends to that found in 
the previous section emerge. In particular, as the aperture 
becomes larger (either in radius or depth), the differences 
between the colour distributions of galaxies in environment 
extremes lessen. Here, the individual properties of galaxies 
are smoothed over due to the large variety of local envi- 
ronments falling within the aperture. For apertures probing 
scales much larger than the typical cluster the distinction 
between environments vanishes. This suggests that environ- 
ment questions relating to galaxy colour (or properties that 
correlate with colour) should avoid fixe d aperture method s 
with large smoothing radii or depth (e.g. lCroton et aLll2005l ). 
Furthermore, the Kolmogorov-Smirnov probabilities indi- 
cate that nearest-neighbour based methods detect stronger 
colour-environment relations than all the apertures tested 
here. 



4.3 Dark Matter Environment 

Dark matter haloes are known to be biased tracers of the 
underlying dark matter distribution, and it is interesting to 
compare how haloes and the smooth background mass field 
correlate with respect to their environment ranking, and how 
this relates to the galaxy distribution. To this end, the sim- 
ulation volume has been divided using a three dimensional 
grid of side-length 2/i~ 1 Mpc, and the neighbourhood dark 
matter density field measured with a Gaussian filter placed 
at the centre of each grid element, smoothed on three dif- 
ferent scales: 2.5, 5.0 and 10/i _1 Mpc (see Section |3~2"|) . We 
compare this to the environment measured directly from 
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Figure 6. Histograms of galaxy colour for the 20 percent most dense (red solid) and 20 percent least dense (blue dashed) galaxies, 
measured using n-th nearest neighbour statistics, defined by the (from left to right) 3rd nearest neighbour, 7th projected nearest 
neighbour, 10 neighbour Bayesian metric and 64 neighbour kernel smoothed (SPH style). The number in the upper right of each panel 
is the Kolmogorov-Smirnov probability that both samples are drawn from the same distribution. 
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Figure 7. Histograms of galaxy colour, as in Figure[6] for the 20 percent most dense (red solid) and 20 percent least dense (blue dashed) 
galaxies, measured using (from left to right) a x Mpc aperture with a velocity cut of ±1000 kms -1 , a 2/i — 1 Mpc aperture with a 
velocity cut of ±6000 kms -1 , a 2/i — 1 Mpc aperture with a velocity cut of ±500 kms -1 and a 5?t — 1 Mpc aperture with a velocity cut of 
±1000 kms -1 . The number in the upper right of each panel is the Kolmogorov-Smirnov probability that both samples are drawn from 
the same distribution. 




Dark Matter Percentage Rank 



Figure 8. The percentage rank of central galaxy environment using an 8 h Mpc spherical aperture plotted against the percentage rank 
of background dark matter environment measured using a smooth Gaussian filter of radius (from left to right) 2.5h -1 Mpc, 5h — 1 Mpc 
and 10/i -1 Mpc. Contours are linearly spaced showing regions of constant galaxy number. 
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central galaxy co unts within a fixed spherical aperture of 
radius 8/i _1 Mpc jCroton et al.ll2005t ). 

Figure [8] shows how the background dark matter den- 
sity, Gaussian smoothed on various scales, correlates with 
the large-scale galaxy density, top-hat smoothed on an 
8/i _1 Mpc scale. The correlation is weakest for the smallest 
Gaussian smoothing scale of 2.5 /i -1 Mpc, becomes tighter 
at a scale of 5/i _1 Mpc, before becoming weaker again at 
10/i _1 Mpc. The point of tightest correlation between dark 
matter and galaxy measured density approximately corre- 
sponds to the same physical scale being probed by each in 
three-dimensional space. At fixed dark matter density the 
scatter in density measured by galaxies is approximately 40 
percent. This indicates the degree of precision with which 
one can probe the smooth background density using galax- 
ies as tracers of the mass distribution. 

We have compared the other environment measures 
used in this paper to the background dark matter density 
but omit the figures for brevity. In short, a similar trend 
to Figure [8] is found for the 64 neighbour smooth density 
environment measure, but with the tightest correlation at 
a radius of 2.5 /i _1 Mpc. For the other neighbour and small 
aperture methods, weak correlations are found when plotted 
against a dark matter density smoothing scale of 2.5 /i -1 Mpc 
but which disappear on larger scales. Environments mea- 
sured in annuli and projected aperture methods that impose 
a photometric-type redshift velocity cut show no correlation 
on any scale due to only the largest clusters dominating the 
depth cut, while the 10 neighbour Bayesian metric and 20 
neighbour smooth density again show a similar correlations 
to the dark matter smoothing scale of 2.5 /i _1 Mpc. 

4.4 Individual Galaxies 

In the previous sections we investigated how different envi- 
ronment parameters correlate with different galaxy proper- 
ties in a statistical sense by considering the whole sample. 
As implied by Figures [3] and [5] when selecting the most 
and least dense environments different methods will poten- 
tially select different galaxy populations. An alternative and 
complementary way to compare the different environment 
methodologies is to consider individual galaxies in the mock 
catalogue and examine how each measure ranks them rela- 
tive to the others. By considering individual galaxies a bet- 
ter understanding of why these galaxies were chosen can 
be obtained. This also highlights the consistency (or lack 
thereof) between different definitions of environment. Below 
we present one example that is representative of the general 
behaviour for high mass haloes. 

The top panel of Figure [9] shows how the different en- 
vironment measures listed in Table [I] compare when one 
focuses on the central galaxy occupying the fourth most 
massive halo in the simulation, with mass io 15 08 h~ 1 MQ. 
The environment measures are separated into three groups 
based on the technique they use: neighbours, aperture and 
annulus. All environment measures place this galaxy within 
the top 10 percent of rank ordered densities in the simu- 
lation volume, with the majority placing it within the top 
1 percent. When considering annuli to define environment, 
the top panel of Figure [9] shows that the further one moves 
from the centre of the halo the lower the rank density mea- 
sured. This simply highlights that the outer regions of a 
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Figure 9. (top) The percentage rank of various density estima- 
tors (see Table [TJ for a single central galaxy living in the fourth 
most massive halo in the simulation. The density methods are 
grouped by increasing neighbour number, increasing aperture and 
increasing inner radius of an annulus. (bottom) The same as the 
top panel, but this time the percentage density ranking of the 
outer most satellite galaxy in the same halo for each method. 



halo tend to be less dense than the core. When considering 
aperture methods there is less of a trend between different 
definitions. However, for a fixed depth, increasing the aper- 
ture size reduces the rank density measured, while for fixed 
aperture size, the density rank appears sensitive to the in- 
clusion of both the halo core (smaller velocity cut) and full 
extent (larger velocity cut). As mentioned in previous sec- 
tions, the larger velocity cut used to represent photometric 
redshift uncertainties has a smaller effect on large clusters as 
the cluster members dominate the galaxies within the depth 
cut. For neighbour-based methods there is a general increase 
in the rank density as the neighbour number increases. This 
is due to the increased neighbour count contributing from 
within the galaxy halo. Specifically, as the number of neigh- 
bours increases galaxies in smaller haloes are demoted down 
the rank list, and so the galaxies in large haloes are pro- 
moted. 

The bottom panel of Figure [9] shows how the different 
density estimators rank the most distant satellite associated 
with the central galaxy of the same 1O 15 O8 M0 halo used in 
the top panel. This is a test of how satellites on the outskirts 
of cluster environments would be classified in each density 
scheme. The range of environment ranking is much larger 
between the different measures, with the apertures and an- 
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nuli mostly finding higher rank densities than the neighbour 
estimates. This again comes back to scale, with neighbour 
methods probing the internal properties of the outer halo, 
and aperture and annulus methods being sensitive to the 
larger structure of the halo and its surrounds. Additionally, 
the trend of increasing rank density with increasing neigh- 
bour number is again seen as the neighbour count reaches 
deeper into the halo core from the boundary. 



5 DISCUSSION AND SUMMARY 

The phrase "galaxy environment" is a very general concept 
that has been used in the literature in a variety of ways. Its 
definition - what it measures and how it is measured - can 
vary from author to author. This creates uncertainty when 
trying to compare results for environmental trends. In prac- 
tise, galaxy environment is quantified in one of two ways: 
by the distance to the n-th nearest neighbour or by using 
a fixed aperture to probe the surrounds. Over the course of 
time these two methods have evolved in the literature. How- 
ever, both methods and their variants provide a measure of 
the density field surrounding a galaxy and hence can be used 
to answer specific environment-related questions. 

To fairly compare many different environment measures 
one would ideally like to use a common galaxy catalogue as a 
starting point. This was achieved in the present work by ap- 
plying a halo occupation distribution model to the z = out- 
put of the Millennium dark matter simulation. Our model 
is designed to accurately reproduce the luminosity, colour 
and spatial distribution of galaxies in the Sloan Digital Sky 
Survey. The resulting data cube was also used to generate a 
mock light cone so that the environment measures could be 
applied in a more realistic geometry. 

The major conclusions of this paper are as follows: 

• Galaxy environment versus halo mass: Compar- 
ing neighbour and aperture based environment measures to 
the dark matter halo mass of a galaxy reveals how they 
measure different aspects of the halo. In particular, near- 
est neighbour methods that use a small enough neighbour 
number best probe the internal properties of the halo. For 
haloes that contain fewer galaxies than the neighbour num- 
ber, the inter-halo separation dominates the calculation and 
galaxy-environment correlations tend to wash out. In con- 
trast, aperture measures tend to better probe the halo as a 
whole and so lead to larger density values corresponding 
to larger haloes, which more accurately reflect their larger 
masses. A smaller aperture than those studied here could 
be used to probe cluster environments on a scale similar to 
the nearest neighbour-based methods, but these would be 
unsuitable for the field due to the distance between neigh- 
b ours being too large. This is in agreem ent with the findings 
of lHaas. Schave fc Jeeson-Daniell (|2012T ). 

• Galaxy environment within haloes: Galaxy density 
internal to a halo's boundary is found to be independent of its 
mass when probed using the neighbour method. While galax- 
ies at the edge of a halo are always in less dense environments 
than those at the centre, the galaxy environment at the cen- 
tre of intermediate mass haloes is approximately the same 



as that at the centre of very massive ones0 By fitting the 
number of galaxies for a given halo mass, we find that the 
number of galaxies per unit dark matter mass is constant 
and th is is in agreement with the findings of lPoggianti et all 
(|20ld ). 

• Environmental dependence of galaxy colour: 

When comparing how the different environment measures 
distribute galaxy colour, almost all methods recover the ob- 
served correlation that galaxies are redder in denser environ- 
ments compared to those in less dense environments. This 
relation only disappears for very large apertures, of order 
> 5 /i _1 Mpc. On scales larger than this the most dense and 
least dense galaxies are found to have similar colour distri- 
butions. Here, the aperture is large enough to encompass 
a statistically representative number of different haloes, re- 
sulting in a smoothing out of the colour differences over such 
large volumes. This behaviour is also expected to extend to 
any property that correlates strongly with colour. 

• Galaxy environment versus dark matter envi- 
ronment: On the other hand, very large fixed apertures are 
the most accurate at recovering the large-scale dark matter 
environment. For example, an 8/i _1 Mpc spherical aperture 
used to calculate the galaxy density correlates well with the 
dark matter environment measured using Gaussian smooth- 
ing on 5/i _1 Mpc scales. Similar results are found with high 
number n-th neighbour estimates. The important parame- 
ter here is scale, with larger probed scales better correlated 
than smaller scales. 

• Environment measures using spectroscopic ver- 
sus photometric redshifts: In addition to the environ- 
ment measures themselves, we also explored the general ef- 
fects of photometric and spectroscopic redshift uncertainties 
by varying the velocity cut used to calculate projected en- 
vironment. For a typical photometric redshift uncertainty 
most trends with environment disappear or become signifi- 
cantly weaker. This is caused by the depth of the aperture 
becoming much larger than the objects being probed, and 
the scattering of interlopers which contaminate the density 
probe. This effect decreases for the largest clusters as the 
members dominate the depth cut. We warn that photomet- 
ric redshifts may be unsuitable for measuring certain proper- 
ties when using a range of environmental scales, especially 
at high redshifts. 

• Environments of central and satellite galaxies: 
On a galaxy-by-galaxy basis, most environment methods 
agree on the relative environment rank of central galaxies 
in massive haloes (to within a few percent). There is less 
agreement with the satellite population in clusters, where the 
result becomes more sensitive to the method employed. 

Dark matter haloes are often broadly categorised as re- 
siding in 'field', 'group', and 'cluster' environments based 
solely on their mass. In Figure [T] for example, haloes of 
mass M ~ 10 12 M Q , 1O 13 - 5 M , and 1O 15 M approximately 
correspond to these environment bins, respectively. Many 



3 The concentration and mass of dark matter haloes are anti- 
correlated, and since the number density distribution of ga laxies 
follows that of the dark matter particles l lYang et alj|2005h . the 
central concentration of galaxies should also vary slightly with 
halo mass. In practice, however, the trend is difficult to detect 
observationally. 
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environment analyses use this categorisation, although as 
we have seen, the distinctions between them can often be 
blurred in detail. Some studies also attempt to explicitly 
identify galaxies that are isolated or in groups or rich clus- 
ters, for exa mple using friends-of-friends group-fin ding algo- 
rithms (e.g. iBerlind et aljlioo^ ; lYang et alj|2007f ). Analyses 
using group catalogues are complementary to studies with 
nearest neighbo ur or fixed aperture mea s ures, or with galaxy 
clust e ring (e.g. Weinmann et all 20061; Martinez fc Muriel 



20061; IBlanton fe Berlind 120071; Ivan den Bosch et al.l l2008aT 



Balogh et all 120091 ; ISkibba et all l201lh . Work focused on 



Poggianti et al.l 


2008; iRudnick et al. 


120091; iBamford et all 


2009; Wolf et al. 


20091; iGallazzi et al. 


2009). 



Importantly, the way a galaxy forms and evolves is 
clearly related to its environment. Some galaxy proper- 
ties, such as luminosity, colour, and stellar mass, are di- 
rectly correlated with the large -scale environment through 
the h ost dark matter halo fe.g. IZehavi et al.ll2005l ; iLi et al.l 
2006). Other galaxy properties, such as structure, are 
to some exte nt only indirectly corre l ated with the envi- 
ronment (e.g. Kauffmann et all |2004| ; IBlanton et al. 1 120051 ; 
ICassata et al.l 20071 ). Indeed, a number of authors report 



that, for many aspects of the galaxy population, envi- 
ronme ntal dependence is often w ea k once stellar mass is 
fixed (Ivan den Bosch et alj l2008allbl; IScodeggio et all [20091 ; 



iBolzonella et al.l l20ld ; IVulcani et all 120111 )/ In any case, 
these studies highlight the fact that it is important to care- 
fully determine how a galaxy's environment is characterised, 
and to identify and navigate the potential aspects of the en- 
vironment analysis that may bias the results. 

The key consideration when picking an environment 
measure is the scale that is being probed. The term environ- 
ment is very general but in fact breaks down into two main 
regions, and we argue that the community should agree on 
a standard terminology for clarity and to avoid future con- 
fusion. The first region is the 'local environment' which cor- 
responds to scales internal to a halo. These are best probed 
using nearest neighbour methods, but the value of n is im- 
portant. When n is larger than the number of galaxies likely 
to reside within the halo the usefulness of this environment 
measure can weaken. The second region lies external to the 
halo, the 'large-scale environment'. The large-scale environ- 
ment is best probed using aperture based methods. In gen- 
eral, there is no simple way to probe all environments with 
a single method, and one should consider carefully the best 
tool to answer the questions at hand. 

This paper marks the first in a series exploring the 
meaning and methods of galaxy environment, as measured 
in the current literature. In the present work we have fo- 
cused on using a clean sample of mock galaxies to quantify 
how selected properties of the galaxy population correlate 
with different environment methods, and how these methods 
themselves compare. Future work will include investigating 
the detrimental effects of survey geometry, edges and holes 
(such as those caused by stars) on environment and tech- 
niques that can be applied to successfully overcome them. 
Furthermore, the relationship between galaxy, halo and dark 
matter environment warrants additional exploration, as does 
the redshift dependence of a galaxy's environment (defined 
in various ways), what the different environment methods 



tell us about galaxy evolution, and how these can best be 
applied to the noisy data of the high redshift Universe. 
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